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1. Introduction 

We present an update of our long term investigation of properties of % and K mesons ||, ||] 
with three flavors of improved staggered quarks [Q]. Study of the % — K system is interesting 
because it enables: 1) a sensitive check of algorithms and methods, including the use of the fourth 
root of the determinant for dynamical staggered quarks, by comparing lattice results for f % to the 
well-determined experimental value; 2) a precise extraction of the CKM matrix element V us from 
fx or fx/ fn, competitive with the world-average from alternative methods; 3) a determination 
of the light quark masses and their ratios with high precision; 4) a calculation of several of the 
physical coefficients of the &{p ) chiral Lagrangian, the Gasser-Leutwyler L,-; 5) a test of the 
applicability of staggered chiral perturbation theory (SZPT) ^ ||,[7|] for describing lattice data; and 
6) a determination of the extra, unphysical parameters that enter SZPT. This allows, for example, 
use of heavy-light SZPT ^ for computations of decay constants and form factors in D and B 
systems []9j] without introducing large additional uncertainties. 



2. Gauge Configurations and Chiral Fits 

The MILC Collaboration has generated an extensive set of gauge configurations with three 
flavors of improved staggered quarks: m' u = m' d = m' , and m' s . (Primes indicate simulation values; 
corresponding masses without primes are the physical values). In this work we report results for 
ensembles with lattices spacings a rs 0. 15, 0. 12, 0.09 and 0.06 fm, which are referred to as coarser, 
coarse, fine and super-fine lattices, respectively. Coarser, coarse and fine lattices are available with 
a range of light quark masses, the lowest pion mass being m n R* 240 MeV. A run on the super-fine 
lattices with m n rj 430 MeV is half-finished, and lighter-mass runs are being started. The simulation 
strange quark masses m' s are in the range 0.70m s <m' s <l.2m s . The physical volumes of the lattices 
range from rs (2.4 fm) 3 to rj (3.4 fm) 3 , with the large volumes being used for ensembles with the 
lightest quark masses. The lattice spacing is kept fixed within each ensemble (coarser, coarse, fine 
and super-fine) as the light quark mass is varied, using the length r\ [JlO, 11] from the static quark 



potential to set the scale. The absolute scale is set from the T 2S- IS splitting determined by the 



HPQCD Collaboration [ 12, 13] on most of our lattices. From their results we find r\ = 0.318(7) fm. 

For Goldstone masses and decay constants, we have extensive partially quenched data, typi- 
cally all combinations of 8 or 9 valence quark masses between 0.1 m^. and m' s . Goldstone quantities 
have the smallest statistical errors, so we concentrate on them. We fit decay constants and masses 
together including all correlations, and we fit different lattices spacings together. The confidence 
level of the joint fit is 0.99. Statistical errors are very small: 0.1% to 0.8% for squared masses, and 
0.1% to 1.4% for decay constants. 

In order to obtain good fits to SZPT forms, we need to place upper limits on quark masses. 
We consider two data sets: 

Subset I: In this subset the valence Goldstone pion masses are restricted to be < 350 MeV, and 
the other pion masses to be < 550 MeV for the coarse ensemble, <460 MeV for the fine ensem- 
ble, and < 400 MeV for the super-fine ensemble. The largest sea quark masses and the coarser 
ensemble are excluded, leaving a total of 122 data points. We expect errors of NLO SZPT to be 
of order 0.8%, implying that NNLO terms are needed. NNLO SZPT logs are unknown; however, 



2 



Update on the physics of light pseudoscalar mesons 



R. Sugar 



for higher masses, where NNLO terms are important, such logs should be smoothly varying and 
well approximated by NNLO analytic terms. The NNLO fit has 20 unconstrained parameters. An 
additional 6 tightly constrained parameters allow for variation of physical LO and NLO parameters 
with the lattice spacing (~ asa 2 AQ CD f» 2%), giving a total of 26 parameters. This subset is used 
to determine the L,, and the systematic errors on other quantities. 

Subset II: In this subset the valence Goldstone pion masses are restricted to be < 570 MeV, and the 
other pion masses to be < 780 MeV for the coarser ensemble, < 725 MeV for the coarse ensemble, 
< 630 MeV for the fine ensemble, and < 590 MeV for the super-fine ensemble. All sea quark 
masses and the coarser ensemble are included, leaving 978 data points. Even NNLO fits break 
down for this subset. We therefore add 18 NNNLO analytic terms and 10 tightly constrained 
variations of NNLO parameters with a. We fix (within errors) LO and NLO terms from the fit to 
subset I. We are then left with 28 unconstrained parameters and 26 tightly constrained ones. This 
data set is not chiral, but it is used to interpolate around the physical value of m s . The central values 
of decay constants and quark masses are determined from it. 

We emphasize that the necessity for using high order SZPT fits with large numbers of parame- 
ters arises from the very small statistical errors in our data. If we did not care about the confidence 
level, we could perform a NLO SZPT fit with only 12 parameters. It changes f n by 4%, by 
1%, and m s by 3% on Subset I; however, it has % 2 /d.o.f. = 9.5 for 110 d.o.f.. Similarly, the 
NLO SZPT fit with taste- violating parameters given as input and the lattice-spacing dependence of 
physical parameters set to zero has 6 parameters. It changes f n by 2%, fx by 1% and m s by 0.5% 
for Subset I; but it has x 2 /d.o.f . = 40.5 for 110 d.o.f.. Furthermore, such fits can change the low 
energy constants by 3 or 4 a. It is important to note that if the physics is not correct, one will not 
be able to fit the data even with ~40 parameters. This was illustrated by a series of fits in Ref. [ffl]. 
There, comparable fits to the continuum XVT form had 36 parameters, but j 2 /d.o.f . = 8.8 for 
204 d.o.f., for a confidence level of 10~ 25() . Fits with all chiral logs and finite volume corrections 
omitted from the fit function (analytic function only) were poor, providing good evidence for chiral 
logs. Finally, separate linear fits for nr % or f n as a function of the average valence quark mass were 
tried. For m| there were 6 parameters and % 2 /d.o.f . « 20 for 234 d.o.f., while for f n there were 
10 parameters, and % 2 /d.o.f . « 25 for 230 d.o.f.. 

3. Taste Symmetry and Chiral Logs 

Violations of taste symmetry for staggered quarks are most apparent in the pion sector. Fig- 
ure [I] shows the splitting between the square of Goldstone pion mass and the square of the masses 
of pions with other tastes for the coarse, fine and super-fine ensembles with m' = 0.4m^,. The split- 
tings are plotted as a function of a 2 a 2 , the variable in which they are expected to be linear at small 
lattice spacings. Figure || shows the square of the pion masses as a function of the light quark mass. 
The degeneracies predicted by Lee and Sharpe in Ref. [Q] are clearly visible. These two figures 
provide strong evidence that taste symmetry violations are well understood from SZPT, and that 
they decrease as expected with decreasing lattice spacing. 

The data is accurate enough and at small enough quark masses so that the effects of chiral 
logs are evident. This is illustrated in Fig. || where the square of the pion mass divided by the sum 
of the valence quark masses m x + m y is plotted as a function of m x + m y in units of r\ . The mass 
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Figure 1: Taste splitting among the pions for lat- 
tice spacings 0.12, 0.09 and 0.06 fm with m 1 = 
0Am' s . 



Figure 2: Squared pion masses as a function of 
the light quark mass on the coarse lattices. 



renormalization constant, relative to that of the fine lattices, has been included so that data from all 
lattice spacings can be presented on the same plot. Lines through the data points come from the 
SZPT fit to the entire data set for decay constants and masses discussed above. The upward slope 
in these lines at small valence quark masses is caused by partially quenched chiral logs, which are 
of the form {m 2 sea G + A 1 fl ,2 a|)log(m^ a/ G + A,a 2 a|), where m va i c and m sea Q are the valence and 
sea Goldstone pion masses, respectively, and A, parametrizes the splitting of the taste-/ pion. The 
effect is especially pronounced in the super-fine data because, not only are the taste- violating mass 
splittings rather small, but the light sea quark mass is rather large, m 1 = 0.4m' s . The red line is the 
fit function in "full continuum QCD" (valence and sea quark masses set equal, and extrapolation 
of the parameters to the continuum limit). It is much smoother because it does not have partially 
quenched chiral logs. 



4. Preliminary Results 

Our previously published work [[jj [2], ||] was based on analysis of the coarse and fine ensembles. 
Here we present preliminary results obtained by including the coarser and super-fine ensembles, 
and additional data from the m' = 0.1 m' s run in the fine ensemble. 

Figure ^ shows fit for the leptonic decay constant of the pion as a function of the sum of the 
valence quark masses, in units of the scale r\ . Lines through the data points come from the same 
SZPT fit as in Fig. [$]. The red lines represent the fit function in "full continuum QCD" with the 
strange sea quark mass fixed to its physical value. The red plus shows the extrapolated value of f % 
to the physical point, in agreement with experiment (black burst) to within systematic errors (blue 
bar). 

Preliminary numerical results for f % and fx obtained from our most recent data are 
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Figure 3: The square of the pion mass divided 
by the sum of the valence quark masses as a 
function of the sum of the valence quark masses 
in units of r\ . 



Figure 4: The pion decay constant as a function 
of the sum of the valence quark masses in units 
of r\. 



A = 128.6 ±0.4 ±3.0 MeV [129.5 ±0.9 ±3.5 MeV] (4.1) 
f K = 155.3±0.4±3.1MeV [156. 6± 1.0±3.6MeV] (4.2) 
/k/A = 1.208(2)(1 14 7 ) [1.210(4)(13)]. (4.3) 

Here the numbers on the left are the new values, and those on the right in square brackets 
are from Refs. [JTJ, @, g. In each case the first error is statistical, and the second systematic. The 
result for f n obtained from the experimental rate for % — > /iv coupled with the value of V u d from 
super allowed nuclear beta decay is f n = 130.7 ±0.1 ±0.4 MeV. The agreement of our result 
with the experimental one provides important evidence that we do understand and can control our 



errors. Marciano has pointed out Jl4| ] that lattice results for fx/ fn can be combined with the 



experimentally determined rate for K^jxv and V lu j to calculate V us . We find 



| V m | =0.2223 (±?6) [0.2219(26)]. (4.4) 

Again, our latest result is on the left with the previous one in square brackets on the right. The 
Particle Data Group (2006) gives V us = 0.2257(21) [JT3J] from the K^Tijiv experimental rate and 
non-lattice theory. Lattice errors continue to dominate experimental ones in our determination of 
Vus, and they will be reduced as additional super-fine lattices become available. 

The up, down and strange quark masses can be determined from the masses of the 71 and K 
mesons using the SZPT fits. To do so we must distinguish between experimental masses, QCD 
masses in which electromagnetism has been turned off, and those in which both electromagnetism 
and isospin violations are turned off. The last of these, which we denote by m% and mg, are used 
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to determine the physical values of the bare quark masses m and m s from the SZPT fits. They are 
related to the QCD masses by 



in 



in, 



in 



(4.5) 
(4.6) 



To relate the QCD masses to the experimental ones requires some continuum input [16]. De- 
tails can be found in Ref. [JTJ] _ Once m and m s have been determined, m u can be obtained from 
m@+ D . Preliminary results from our current data set are: 



m,. 



■ 90(0) (5) (4) (0) MeV 
3.3(0) (2) (2) (0) MeV 
m s /m = 27.2(0) (4) (0)(0) 
f^-=2.0(0)(l)(l)(l)MeV 
m MS = 4 6 ( )( 2 )(2)(l) MeV 
m M /m d = 0.42(0)(l)(0)(4) 



m 



MS 



[76(0) (3) (7) (0) MeV] 
[2.8(0)(1)(3)(0) MeV] 
[27.4(1) (4) (0)(1)] 
[1.7(0)(1)(2)(2) MeV] 
[3.9(0)(1)(4)(2) MeV] 
[0.43(0)(1)(0)(8)] 



(4.7) 
(4.8) 
(4.9) 
(4.10) 
(4.11) 
(4.12) 



Again, new results are on the left, and early ones in this case from Refs. [fiL 13, 17] are in 
square brackets on the right. Errors are from statistics, simulation systematics, perturbation theory 
and electromagnetic effects. The MS scale is /I =2GeV. The main difference between the new 
and old results comes from the perturbation theory calculation of the renormalization constants 
needed to match the lattice masses to the MS masses. Here we use the new two-loop results of 



Mason, Trottier and Horgan [18], whereas the old results were based on a one-loop calculation. A 
non-perturbative mass renormalization calculation is in progress. Note that our result for m u jm,i 
rules out a massless m u at the 10 a level. 

Finally, we have determined a number of the low energy constants L, at the scale of . Our 
preliminary results, in units of 10~ 3 , are: 



2L 6 -L 4 = 0.5(1)(2) [0.5(2)(4)] 

2L 8 -L 5 = -0.1(1)(1) [-0.2(1)(2)] 

L 4 = 0.1(2)(2) [0.2(3) (3)] 

Ls=2.0(3)(2) [1.9(3)(3)]. 



(4.13) 
(4.14) 
(4.15) 
(4.16) 



Again, the new values are on the left, and the older ones from Ref. ^ are in square brackets 



on the right. These results are consistent with "conventional" ones summarized in Ref. L5 = 
2.2(5), Ld = 0.0(3), and L4 = 0.0(5). Our value for 2Lg — L5 is far from the range —3.4 < 2Lg — 
L5 < —1.8 that would allow m u = 0, which is consistent with, but not independent of, our direct 
determination of m u . Our old results for the low energy constants already play a significant role in 
setting constraints on the n — n scattering lengths and Aq [20]. It will be interesting to see the 
impact of the new ones. 
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